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The problem of two magnetic impurities in a normal metal exposes the two opposite tendencies 
in the formation of a singlet ground state, driven respectively by the single-ion Kondo effect with 
conduction electrons to screen impurity spins or the Ruderman-Kittel-Kasuya-Yosida interaction 
between the two impurities to directly form impurity spin singlet. However, whether the competition 
between these two tendencies can lead to a quantum critical point has been debated over more 
than two decades. Here, we study this problem by applying the newly proposed natural orbitals 
renormalization group method to a lattice version of the two-impurity Kondo model with a direct 
exchange K between the two impurity spins. The method allows for unbiased accessing the ground 
state wave functions and low-lying excitations for sufficiently large system sizes. We demonstrate 
the existence of a quantum critical point, characterized by the power-law divergence of impurity 
staggered susceptibility with critical exponent 7 = 0.60(1), on the antiferromagnetic side of K when 
the interimpurity distance R is even lattice spacing, while a crossover behavior is recovered when R 
is odd lattice spacing. These results have ultimately resolved the long-standing discrepancy between 
the numerical renormalization group and quantum Monte Carlo studies, confirming a link of this 
two-impurity Kondo critical point to a hidden particle-hole symmetry predicted by the local Fermi 
liquid theory. 

PACS numbers: 75.10.Hf,71.15.Dx,75.20.Hr,75.40.Cx 


I. INTRODUCTION 

It is well-recognized that the competition between 
the single-ion Kondo effect and the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction, the two inevitable 
forces in any Kondo systems with more than one local 
impurity magnetic moment, plays a crucial role in cor¬ 
related systems ranging from dilute magnetic alloys to 
heavy fermion compounds 0 !- The indirect RKKY 
interaction, namely the interimpurity interaction medi¬ 
ated by conduction electrons via a short-range (on-site) 
Kondo coupling, oscillates and decays with the inter¬ 
impurity distance R and Fermi momentum 2fcp 00. 
When the RKKY interaction grows toward the strong 
antiferromagnetic limit, the quantum many-body ground 
state will evolve from a collective Kondo singlet state [6!] 
into the interimpurity singlet state locked by the RKKY 
interaction. However, whether a distinct separation or 
a quantum critical point exists between the two singlet 
ground states has remained elusive. This theoretical is¬ 
sue has a fundamental importance as it closely correlates 
with the critical divergence or scaling behavior of several 
physical quantities, or emergent energy scales in realistic 

materials |3-[9( ■ 

It is remarkable that even for the simplest case with 
only two local impurities, i.e., the two-impurity Kondo 
model (TIKM)0, the evidences for a quantum critical 
point separating the two distinct singlet states do not 
converge. On the one hand, the early numerical renor¬ 
malization group (NRG) studies flll4l3l ] revealed an un¬ 
stable fixed point characterized by diverging impurity 


staggered susceptibility and coefficient of specific heat 
at a finite ratio of the RKKY interaction to the Kondo 
temperature (« 2.2). On the other hand, the subsequent 
well-controlled quantum Monte Carlo (QMC) studies did 
not find such a divergence, instead, a crossover behavior 
in the corresponding quantities at very low but still finite 
temperatures was observed 00. Since then, although 
the effective field theory analysis suggested the occur¬ 
rence of critical points in several variants of the TIKM 
0-0, the strong debate in the TIKM studies, espe¬ 
cially in numerical studies, has still remained until today 

nil. 

Almost all the previous investigations, except for the 
QMC studies 0,0, rely on a decomposition of the con¬ 
duction electrons into odd/even channels with respect to 
the impurity center. Under such a decomposition, the 
original two-impurity problem is mapped effectively onto 
two-channel/two-impurity problem defined on the mo¬ 
mentum or energy space, resulting in various impurity 
couplings which are energy-dependent in general. The 
obtained two-channel/two-impurity model, capturing the 
low-energy properties of the original problem, provides a 
base for the NRG and effective field theory studies. In 
particular, the early NRG calculations 0, 0 assume 
“energy-independent” coupling constants in odd/even 
channels. The QMC simulations 00 , on the other 
hand, suffer from finite temperatures. Other NRG [23- 
0 or DMRG0 studies, where no evidence was found 
for a quantum critical point after taking into account the 
“energy-dependent” coupling strengths or starting from 
a real-space two-impurity lattice model, indeed show the 
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importance of lattice geometry details. 

More seriously, these numerical results seem to 
strongly contradict a general phase shift argument based 
on the local Fermi liquid theory [2(1 [23], which states that 
a phase transition must exist between the two stable fixed 
points if a TIKM preserves particle-hole (PH) symmetry 
[281 f29j | . As a matter of fact, there are two types of 
PH symmetries associated with a standard TIKM, cor¬ 
responding to the cases with the interimpurity distance 
R being even or odd, or in the lattice case, the impurity 
mirror center being on site or bond, respectively. It is the 
first type of PH symmetry, namely the distance R being 
even or the mirror center being on site, that can guaran¬ 
tee a phase transition mu. Although the phase shift 
argument comes from the odd/even decomposition with 
a number of simplifications including the spherical plane 
wave and linear dispersion approximations, its validity 
should not depend on ignoring energy dependence of cou¬ 
pling constants and other lattice details. Otherwise, if 
there is no any phase transition in this case, it would in¬ 
dicate a breakdown of the local Fermi liquid picture of 
the single-impurity Kondo problem. Here we would like 
to emphasize that such a phase transition could be either 
a first-order transition or a quantum critical point. In the 
case of only the second type of PH symmetry preserved, 
there is no such a guarantee, in other words, either a 
phase transition or a crossover takes place between the 
two stable fixed points. 

In our point of view, the most straightforward ap¬ 
proach to challenge or confirm this argument and fur¬ 
ther clarify the discrepancies among various numerical 
studies is to directly solve the ground state of a stan¬ 
dard TKIM for sufficiently large systems without using 
decomposition. In fact, Affleck et al. [29j] already outlined 
three conditions for such a decisive numerical study: (1) 
the studied model preserves the first type of PH symme¬ 
try; (2) a model parameter is varied to pass the interim¬ 
purity singlet to the Kondo singlet; (3) sufficiently low 
temperatures are accessible and other model parameters 
are fine-tuned. 

In this paper, we reexamine whether or not there is a 
quantum critical point in a standard TIKM by applying 
a newly developed numerical method, i.e., the natural or¬ 
bitals renormalization group (NORG) HI} . Different from 
the conventional NRG, the NORG keeps faithfully all the 
lattice geometry details and does not need a decomposi¬ 
tion into odd/even channels and mapping onto the mo¬ 
mentum or energy space. In particular, the ground state 
wave function as well as the relevant low energy excita¬ 
tions in electronic bath with several impurities can be 
calculated accurately in a well-controllable manner for 
a very large system. Therefore, the NORG is particu¬ 
larly pertinent for studying the zero temperature quan¬ 
tum phase transition in the two-impurity Kondo prob¬ 
lem. 

In the following we will present the numerical results 
and discuss the main observations. A brief description of 
the NORG method and an illustration of the phase shift 


FIG. 1: (Color online) Schematic view of a standard two- 
impurity Kondo model with impurity separation R = 3. 

argument and other relevant numerical results are given 
in the appendixes. 

II. MODEL AND RESULTS 

Our studied TIKM is a standard one-dimensional lat¬ 
tice model 52], as schematically shown in Fig. [TJ which 
is described by the following Hamiltonian, 

H = H 0 + F^Kondo + Hrkky; 

Ho ^ y~)[0/r0-t-ifT T h.c.], 

HKo„do = J £ Sj ■ s (Rj), (1) 

-Bricky = KSi ■ S 2 , 

where Hq describes the non-interacting conduction band 
with Ci a being the annihilation operator of a conduction 
electron located at the i-th site with spin component a 
and t being the nearest-neighbor hopping integral, J is 
the short-range (on-site) Kondo coupling between each 
impurity spin S ; and the spin s(Rj) of a conduction 
electron passing by the j’-th impurity site, and K is 
the direct exchange interaction between the two impu¬ 
rity spins. Such a TIKM can be realized in nanoscale 
devices where the observed Kondo signature varies with 
tunable RKKY interaction [32435} . Notice that in real¬ 
istic materials the direct interaction K, which is finite 
even at vanishing Kondo coupling, can be mediated as a 
superexchange via othe_r_filled orbitals as in some layered 
/-electron compounds [361 ■ 

As shown in Fig. [I] the two impurities are located 
R = |R 2 — Ri | lattice spacings apart. The model exhibits 
the lattice inversion symmetry with respect to the mirror 
center of the two impurities. As illustrated by Affleck et 
al.jlsj, a t half filling of the conduction band the model 
exhibits the first (or second) type of PH symmetry when 
R is even (or odd). Throughout the present study we 
fix t = 1/2 and take half filling of the conduction band, 
while the periodic boundary condition is imposed and the 
length of the chain is denoted by L. When K = 0, the 
model becomes the one studied by Fye et al. using QMC 

method [ME!. _ 

Motivated by the suggestion of Affleck et al.[29;], we 
mainly consider the two representative cases, i.e., the 
interimpurity distance being odd (R = 1,3) or even 
(R = 2,4), respectively. We note that the case with 
even R was not thoroughly considered in the previous 
QMC calculations [3, [n|. However, the required first 
type of PH symmetry is realized only in this situation. 
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We also notice that actually in most of the previous nu¬ 
merical studies iiHi only indirect RKKY 
interaction is considered, corresponding to AT = 0. 

In the present work, the eigenvalues and wave functions 
of the ground state and low-lying states of Hamiltonian 
m were first directly solved by the NORG method for 
sufficiently large systems. The relevant physical quan¬ 
tities were then calculated. Here the proper physical 
quantities to characterize the ground state structure of 
the system and describe its response to external field are 
respectively the interimpurity spin-spin correlation func¬ 
tion (Si ■ S 2 ) ((•••) denotes ground state expectation 
value) and impurity staggered susceptibility Xs that is 
defined as Xs = / 0 °° dr ([Sf (r) — Sf ( r )] [^i ~ Sf]) with 
Sj(r) = e rH Sfe~ Tli . To monitor the change of physical 
property of the system upon varying couplings K and 
J, we thus calculated (Si • S 2 ) and Xs as functions of K 
and J. The corresponding quantum phases were achieved 
by extrapolating (Si ■ S 2 ) and Xs to the thermodynamic 
limit, namely 1/L —> 0, as described in Appendix C for 
details. 

As a benchmark test, we carried out the calculations 
for a set of parameters with I\ = 0, R = 1, J = 0.728 
and 0.8, and with varying lattice sizes L, respectively. 
Figure [2] shows the calculated impurity staggered sus¬ 
ceptibilities along with finite size extrapolations to the 
thermodynamic limit. As we see, the results in the ther¬ 
modynamic limit fall within the error bars of those cal¬ 
culated by the early QMC [HJ. This also confirms the 
early QMC results. We remind that the case of K ^ 0 
was not studied by QMC calculations since this would in¬ 
troduce technical difficulties in QMC. By contrast, there 
is no extra difficulty by using the NORG approach, as 
shown below. 



1/L 

FIG. 2: (Color online) Finite size extrapolation of impurity 
staggered susceptibility \s at K = 0 and at J = 0.728 and 
0.8 respectively for R = 1. Here L = 22, 24, 30, ..., 510, and 
512, respectively. The endpoint at 1/L = 0 of each curve is 
determined by a quadratic polynomial fit of the leftmost four 
points of that curve calculated by NORG. Here the NORG 
results are consistent well with those of Fye ’s QMC [IH|. 



FIG. 3: (Color online) Impurity staggered susceptibility 
and interimpurity spin-spin correlation (Si • S 2 ) (inset) as 
functions of K for J = 1 and R = 2 with different lengths L. 
It is similar to the case of J = 1 and R = 4 (not shown). For 
L = 4n+2 represented by solid lines [37j], Xs shows a divergent 
peak at K = K° dd (L) which moves to the right and finally 
converges to K c = K C (L = 00 ) = 0.1555 as L — > oo[3(|. And 
for L = 4n represented by dotted lines, Xs shows a peak (with 
discontinuity) at K = A'® ven (L) which moves to the left and 
finally converges to the same K c (determined in the L = 4n+2 
case) as L — 00 and meanwhile the peak height increases and 
finally diverges. (Si • S 2 ) jumps where Xs peaks. The dis¬ 
continuity of (Si ■ S 2 ) indicates a quantum phase transition, 
which is further characterized as a quantum critical point by 
the divergence of y s . 


We first study the cases of interimpurity distance 
R = 2 and 4 with Kondo coupling J = 1 respectively, 
in which there is the first type of PH symmetry. Figure [3] 
shows the calculated impurity staggered susceptibility Xs 
and interimpurity spin-spin correlation (Si • S 2 ) at zero 
temperature as functions of K for various lattice sizes. 
Physically a very large negative AT, namely a large ferro¬ 
magnetic coupling, tends to lock the two impurity spins 
into a triplet, characterized by (Si -S 2 ) being positive and 
close to the triplet value j. This has been shown in the 
inset of Fig. [3] The resulting triplet is then screened by 
the conduction electrons through the two-stage Kondo ef¬ 
fect processpjQl]. In particular, we see that this tendency 
extends to AT = 0. This means that the indirect RKKY 
interaction, which is due to the second or higher order 
corrections to the expectation value of (Si • S 2 ) between 
two impurity spins (when AT = 0), is ferromagnetic at 
half filling. Thus the ground state for I\ = 0 is always 
in the Kondo regime. In order to reach the border of the 
interimpurity singlet state regime, we tuned the strength 
of K in a wider range from ferromagnetic to antiferro¬ 
magnetic. 

For a large positive K , namely a large antiferromag¬ 
netic coupling, the ground state shall escape from the 
Kondo screening phase, entering the interimpurity spin 
singlet phase. This is clearly seen from the interimpu- 
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rity spin-spin correlation (Si • S 2 ) (the inset of Fig. ©, 
which is negative and approaching the singlet value — | 
on the large positive K side. In this case, there is the 
development of interimpurity singlet driven by the an- 
tiferromagnetic correlation. In the intermediate regime 
of K ~ 0.15, however, we find a rather sharp peak in 
the impurity staggered susceptibility. With increasing 
lattice size, the peak goes diverging and the position of 
the peak converges to K c , indicating that the divergence 
of impurity staggered susceptibility is an intrinsic prop¬ 
erty free of the finite size effect. Furthermore, when K is 
around the location of the impurity staggered susceptibil¬ 
ity peak, the corresponding interimpurity spin-spin cor¬ 
relation (Si • S 2 ) exhibits a concomitant jump across the 
onset of antiferromagnetic correlation signed by (Si • S 2 ) 
beginning to be negative. This feature indicates a quan¬ 
tum phase transition taking place at K = K c from the 
Kondo singlet regime to the interimpurity singlet regime, 
which is further characterized as a quantum critical point 
by the divergence of the impurity staggered susceptibil¬ 
ity. 

The location of the critical point I\ c can be precisely 
determined by fitting the location of the divergent peak 
in the large L limit, as illustrated in Appendix C. Figure 
U shows the critical point K c as a function of Kondo cou¬ 
pling J, which defines a phase boundary separating the 
Kondo singlet phase (K < K c ) from the interimpurity 
singlet phase (K > K c ). As we see, the phase diagrams 
for interimpurity distance R = 2 and R = 4 are basically 
the same, while the function K C (J) shows a quadratic 
behavior at the small J limit and linear behavior at the 
large J limit. 

To quantify the quantum criticality, we elaborately cal- 



FIG. 4: (Color online) Phase diagram of direct exchange inter¬ 
action K versus Kondo coupling J: calculated critical point 
K c as a function of J. This defines a phase boundary sepa¬ 
rating the Kondo singlet phase (small K) from the interim¬ 
purity singlet phase (large I\) at zero temperature. I\ C (J) is 
quadratic at the small J limit (lower inset) and linear at the 
large J limit (upper inset). 
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FIG. 5: (Color online) Finite-size data collapse for the case of 
J — 1 and R = 2 with different lengths L. This gives rise to 
\ s ~ \K — K c \ y with 7 = 0.60 ± 0.01. Lo = —9.2 and v = 1 
are determined by fitting data with formula \K C (L) — K c \ ~ 
(L — Lo)- 1 /" 0. 


culated the impurity staggered susceptibility \ S (K) as a 
function of K in a range centered at the critical point 
K c with different system sizes, reported in Fig. [5] Us¬ 
ing a standard numerical techniaue(38l|. we can derive 
that Xs{K) diverges towards K c in power law rather 
than logarithmically in the thermodynamic limit, namely 
Xs(K) ~ \K-Kc \~ 1 with 7 = 0.60±0.01. Hence the crit¬ 
ical divergence is faster than the logarithmic one but still 
slower than \K — K c | , indicative of a continuous phase 

transition in computational studies [38|. This exponent 7 
has not been reported in literatures to our knowledge. 

As a comparison, we likewise study the cases of in¬ 
terimpurity distance R = 1 and 3 with Kondo coupling 
J = 1 respectively. There is now the second type of PH 
symmetry, which however does not guarantee any phase 
transition. In Fig. [ 6 ] we report the calculated impurity 
staggered susceptibility \s and interimpurity spin-spin 
correlation (Si • S 2 ) as functions of the direct exchange 
interaction I\. Similar to the cases of R = 2 and 4, for 
a large negative (positive) K, the (Si • S 2 ) approaches 
the triplet (singlet) value, namely - (— |) (the inset of 
Fig©, corresponding to the Kondo (interimpurity) sin¬ 
glet regime. However, there is no any jump or disconti¬ 
nuity for (Si • S 2 ) upon K varying. Actually, as shown in 
the inset of Fig. © the change of (Si • S 2 ) from the Kondo 
singlet regime to the interimpurity singlet regime is very 
smooth across the middle value —7 between the singlet 
and triplet values, indicating co-existence of the Kondo 
and interimpurity singlets in this intermediate regime, in 
other words, a crossover taking place. This can be fur¬ 
ther confirmed by examining the behavior of impurity 
staggered susceptibility Xs as a function of I\ as follows. 
In addition, the (Si • S 2 ), being negative at K = 0, shows 
that the indirect RKKY interaction is antiferromagnetic 
at K = 0, contrary to the cases of R = 2 and 4. 

From Fig. [S]we see that the Xs(K) curves for L = 4n 
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FIG. 6: (Color online) Impurity staggered susceptibility \s 
and interimpurity spin-spin correlation (Si ■ S 2 ) (inset) as 
functions of K for R = 1. J = 1 is fixed with different 
lengths L. For L = 4n + 2, looks diverging at some value 
of K = K c . But the diverging tendency is suppressed and K c 
moves to —00 as L — > 00 [30|,l37|. Therefore, the \s curves for 
L = 4n + 2 and L = 4n both converge to the same one in the 
thermodynamic limit n — > 00 . The case of R = 3 is similar 
(not shown). 


converge to a smooth and non-divergent one as L —>• 00 
M- In contrast, each Xs{K) curve for finite L = 4n + 2 
seemingly has a divergent peak. However, further cal¬ 
culations [30| show that the value of K at which \ S {K) 
diverges for L = An+ 2 moves to the negative infinity as 
L —» 00 . It turns out that eventually the y s (/d) curves for 
both L = 4n+2 and L = An converge to the same smooth 
and non-divergent one in the thermodynamic limit. This 
indicates the absence of a quantum critical point. This 
conclusion is in agreement with the previous QMC cal¬ 
culations for the case of R = 1 (or R = 3). 


III. DISCUSSION AND OUTLOOK 

The studied TIKM offers a new ingredient, namely 
the interimpurity exchange K , which competes with the 
single-ion Kondo effect of individual impurities. At zero 
temperature, two different regimes of the ground state, 
i.e., the Kondo singlet and the inter-impurity singlet, are 
respectively realized when K is varied from the large fer¬ 
romagnetic to the large antiferromagnetic. Crossover be¬ 
tween these two regimes is expected to be a generic fea¬ 
ture since the relevant degrees of freedom due to quan¬ 
tum impurities are finite. Thus it comes as a surprise 
when an unstable non-Fermi liquid fixed point or a quan¬ 
tum critical point separating the two regimes was firstly 
evidenced by the NRG studies [0, Qj]]. The existence 
of such a critical point in the TIKM has been then de¬ 
bated by various approaches even in the presence of PH 
symmetry (13. Qjl [22| - l25j . 


Using the NORG calculations we have demonstrated 
that there is indeed a quantum critical point for a TIKM 
with the first type of PH symmetry. The criticality is 
characterized by the power-law divergence of the impu¬ 
rity staggered susceptibility y s , concomitantly accompa¬ 
nied by a jump of the interimpurity spin-spin correlation 
(Si ■ S 2 ), and thus belongs to a genuine impurity quan¬ 
tum phase transition where the impurity contribution to 
the ground state energy becomes singular [39l-l4l| . More¬ 
over, as presented in above sections, we have also shown 
that to resolve this longstanding issue depends on not 
only the required PH symmetry but also the tunable di¬ 
rect RKKY interaction. Meanwhile, the position of crit¬ 
ical point is closely related to the lattice details. The 
newly developed NORG method, which keeps lattice de¬ 
tails and is free of decomposition and finite temperatures, 
makes examination on all these aspects possible. 

Here we would like to remind that in the previous NRG 
studies on TIKM (or the corresponding two-impurity An¬ 
derson model) the occurrence of a critical point requires 
vanishing even/odd channel asymmetry [24] or energy- 
independence of coupling constants [22], {23[. However, 
the even/odd channel asymmetry or energy-dependence 
of coupling constants is an inevitable consequence of de¬ 
composition when mapping a lattice TIKM onto a contin¬ 
uous model in momentum or en ergy sp ace. On the other 
hand, the phase shift argument [161121.l29j for a quantum 
phase transition in a TIKM does not reply on even/odd 
channel asymmetry or energy-independence of coupling 
constants. Therefore, our NORG study on a standard 
lattice TIKM directly confirms the long-cherished predic¬ 
tion of the link between the first type of PH symmetry 
and the two-impurity Kondo critical point [29j . 

The present study clearly shows that the previous nu¬ 
merical studies actually do not contradict each other 
for their own sake because the discrepancies are mainly 
due to the fact that they adopted different parameter 
regions or approximations. In particular, the energy- 
dependence of the coupling constants appeared in these 
studies should not be an obstacle for studying the critical 
point provided that the first type of PH symmetry and 
antiferromagnetic interimpurity interaction are realized. 
As the required symmetry occurs at half filling in a lattice 
model with reflection symmetry about a site (or interim¬ 
purity distance is even), the two-impurity Kondo critical 
point should be observed experimentally in realistic ma¬ 
terials with tunable interimpurity exchange interaction. 

Finally, the NORG method developed here can be also 
used to solve the multi-impurity models in two or three 
dimensions [32|, or used as impurity solvers in the dy¬ 
namical mean-field theory. Therefore we expect that the 
NORG method can be used to study other challenging 
problems in correlated electron systems. 
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Appendix A: Method 

For a quantum impurity model, by generalizing quan¬ 
tum renormalization group approach into natural orbitals 
space through iterative orbital rotations, we can realize a 
numerical many-body approach, namely the NORG[3l|. 
with polynomial, no longer exponential, computational 
complexity (0(N* ath )) in the number of electron bath 
sites ( Nb a th )■ It turns out that dozens, or even hundreds, 
of bath sites can be dealt with in practice. In addition, 
the NORG works in a Hilbert space constructed from a 
set of natural orbitals, thus it can work on a quantum im¬ 
purity model with any lattice topological structure. In 
the present work, we have improved the efficiency of the 
NORG method. The computational complexity is fur¬ 
ther reduced from 0(N* ath ) to 0{N$ ath ), in which the 
impurity orbitals are no longer involved into orbital ro¬ 
tations, and hence the transformed Hamiltonians become 
simpler. This enables us to calculate the eigenvalues and 
wave functions of ground states and low-lying excitations 
for larger system sizes. The largest system size we have 
reached is L = 1022 (see Appendix C). Errors of the 
NORG are controllable. In this paper the relative errors 
of all data calculated by the NORG are less than 10~ 3 . 


Appendix B: Phase shift argument and particle-hole 
symmetry 

We reiterate the phase shift argument for a phase 
transition between the collective Kondo screening sin¬ 
glet state and the R.KKY interaction-locked interimpu¬ 
rity singlet state in the presence of the first type of 
particle-hole (PH) symmetry (l7l [29t | . In the field the¬ 
ory treatment of the two-impurity Kondo model (TIKM), 
a decomposition of electron field into even/odd compo¬ 
nents ( ip e ,E, V’o.-e) is introduced (the spin degree of free¬ 
dom is implied), both components are dependent on the 
quasi-particle energy E measured from the Fermi energy 
Ep which is set to be zero. The obtained effective Hamil¬ 


tonian, H = Hq + Hk , takes a general form as follows, 


H 0 = 

1 dEE jt l’l, E tpe,E + V’J.eV’o.-e} . 

(Bl) 

H k = 

1 dEdE' ^giip\ E 5ipe,E' ■ (Si + S 2 )J 

' (B2) 

+ . 

1 dEdE' |ff 2 ipl, E v4>o,E' ■ (Si + S 2 )} 


+ . 

1 dEdE' ^g 3 tpl E aip 0! E' • (Si - S 2 )J 


+ 

f dEdE' ^giipl E aip e ^ E ' • (Si - S 2 ) j 



where gi (i = 1,2, 3,4) are energy-dependent coupling 
constants reflecting the lattice details. The first type of 
PH transformation of the fields implies 

V’e.B -A ^l-E 

o,E -A i’l-E, (B3) 

while the second type of PH symmetry implies 

lpe,E -A 1p\ _ E 

Ipo,E -A 1p\_ E - ( B4 ) 

In the Fermi liquid picture, the zero temperature fixed 
point can be characterized by the phase shifts (S e ,S 0 ) 
for both components at the Fermi energy E = 0. This 
amounts to the following relations between the incoming 
and outgoing operators 

r e % = e 2i5 °iP™ E 

C E = e 2iSo <E- (B5) 

Hence 5 e = S 0 in the first type of PH symmetry, and 
S e + S 0 = 0 in the second type of PH symmetry. In the 
limit of K —» — oo, both channels are in the Kondo screen¬ 
ing phase, indicating 5 ei0 = ±7r/2. While, in the limit of 
K oo, the Kondo effect is completely suppressed by 
interimpurity correlation, leading to S e = 5 0 = 0. There¬ 
fore, by varying — oo < K < oo, a phase transition must 
take place at certain value of K c when the first type of 
PH is maintained, as sketched in Fig. [Tfb). By contrast, 
no phase transition is guaranteed for the second type of 
PH symmetry (Fig. [ZKa)) 

Here we would like to remind that for a lattice model 
with inversion symmetry and specified Kondo coupling 
at half filling, the first or second type of PH symmetry 
will be preserved when the interimpurity separation R 
is even or odd, respectively. In general, a long-range hy¬ 
bridization or Kondo coupling can be introduced in order 
to maintain the required PH symmetry. The observation 
described in the main text indicates that a TIKM with 
the simplest on-site Kondo coupling which preserves the 
required PH symmetry is enough to support an emergent 
quantum critical point. 
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Kondo singlet interimpurity singlet 

FIG. 7: (Color online) Schematic view of the phase shift <5 e ( 0 ) 
in the even (odd) channel of the conduction band at the Fermi 
energy. 



1/L 

FIG. 8: (Color online) (Si ■ S 2 ) for R = 1, J — 0.728, and 
K = 0 at zero temperature. L = 22, 24, 30,..., 1022. The end¬ 
point at 1/Z, = 0 of each curve is determined by a quadratic 
polynomial fit of the leftmost four points of that curve calcu¬ 
lated by the NORG. Our NORG result is consistent well with 
that of Fye ; s QMC0. 


Appendix C: Complementary numerical results 

Impurity spin-spin correlation. Another benchmark 
test for our NORG method is the impurity spin-spin cor¬ 
relation, (Si • S 2 ), in the ground state for R = 1, where 
the QMC result is available. We choose J = 0.728 and 
K = 0, the same parameters as given in Fye’s QMC 
study fl5j. The calculated result is shown in Fig. [8] 
The system size varies as L = 22, 24, 30,..., 1022. The 
endpoint at 1/L = 0 of each curve is determined by a 
quadratic polynomial fit of the leftmost four points of 
that curve calculated by the NORG. Our NORG result 
is consistent well with that of Fye ’s QMC[H. The rela¬ 
tive difference between the two endpoint values is just as 
small as 1.4 x 10 -4 , which highlights the high accuracy 



1/J 


FIG. 9: (Color online) K c /T k versus 1/J. In the NRG 
studies H 2 I. Il3l]. K c /Tk was predicted to be a constant, about 
2 . 2 . 

of the NORG method. 

Kc/Tk when R = 2 and 4. To compare with the early 
NRG result, we calculate the ratio K c /Tk, where Tk is 
the single-ion Kondo temperature defined asjlbj 

T k = D{pJ) 1/2 exp(-l /pj), (Cl) 

where D (= 2 here) is the conduction electron bandwidth 
and p (= 1 /7r here) the conduction electron density of 
states at the Fermi level. K c /Tk is then plotted in Fig. [51 
While the generic behavior is similar for R = 2 and R = 
4, quantitative difference between them appears when 
J < 5. It is seen that in a wide range of J, K c /Tk is close 
to the value predicted by the NRG studies, K c /Tk ~ 

2.2.[H[ll] 

Determination of K c when R is even. For even R the 
TIKM has a quantum critical point at K = I\ c where the 
impurity staggered susceptibility Xs diverges. For finite 
system size L = An + 2, \s diverges at K = K C (L). We 
extract these K C (L ) from \s curves calculated by NORG 
for different values of L and extrapolate it to L = 00 by 
fitting. To present the fitting details we show the typical 
case of J = 1 in Fig. [Tl] The perfect fitting shows that 
| K C (L) — K c | ~ (L — To) -1 / 1 ' with u = 1. For a series of 
values of J, we show the result in Fig. [T2] While J spans 
across three orders of magnitude, K c across seven orders 
of magnitude. This shows the robustness of the NORG 
method. 

Suppression of susceptibility peak when R is odd. The 
impurity staggered susceptibility Xs shows a sharp peak 
in the case of R = 3 (or R = 1) when L = An+2. This fea¬ 
ture seems different to what happens when L = An where 
no sharp peak is detected. In order to clarify whether this 
is an intrinsic property of the ground state in the ther¬ 
modynamic limit L —> 00 , we plot the peak position K c 
by varying the lattice size L for different values of J (Fig. 
Hoi. The result suggests that in the limit L —> 00 , the 
peak position K c of Xs for L = An + 2 moves to negative 
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FIG. 11: (Color online) The value of K = K C (L) at which 
diverges for R = 2, J ®= 1, and L = 4n + 2. The dots are 
calculated by NORG. The line is determined by fitting the 
dots with formula K C (L) = K c — a(L — To) -1 , where K c , a, 
and Lo are fitting parameters. This perfect fitting shows that 
| K C {L) - K c | ~ (L - Lo)- 1 /" with v = 1. 


FIG. 10: (Color online) The value of K — Kc at which \ s 
diverges for R = 3 and L = An + 2. Dots denote the results 
calculated by the NORG. Lines are determined by fitting the 
dots with the formula I\ c = —a(L — c) b , where a, b, and c are 
fitting parameters. The K c behavior for R = 1 is similar (not 
shown). This figure clearly shows that K c decreases to — oo 
as L —» oo. 


infinity and % s ’s for both L = An + 2 and L = An con¬ 
verge to the same one with a non-divergent behavior in 
the thermodynamic limit. Therefore, no phase transition 
takes place in the case of R = 3 (or R = 1). 
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